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EXECUTIVE  SUMMARY 


PROBLEM 

Examinc?three-dimensional  matched  filtering  and  show.how  it  can  be  applied  to 

r  n 

accomplish  target  detection  over  sequences  of  observed  scenes,  which  is  an  important 
extension  of  classical  two-dimensional  image  processingf’Provide  examples  to  demonstrate 
and  validate  these  results.  ,  ,  > .  >  *  4  • 


RESULTS 

J Three-dimensional  filtering  techniques  can  be  easily  modified  to  act  as  a  space-shift, 
time-delay-and-integrate  matched-filtering  technique.  Computer  simulations  verified  the 
validity  of  this  modification.  The  results  demonstrate  signal-to-noise  ratio  gains 
asymptotically  growing  by  the  filter  integration  time  over  those  normally  obtained  by 
two-dimensional  matched  filtering  of  the  individual  frames  in  a  sequence.  The  specific  form 
of  the  output  is  ideally  suited  for  moving-target  indicator  (MTI)  applications,  or  can  be 
used  in  conjunction  with  some  other  target-speed  and  detection-localizing  algorithm. 
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1.0  INTRODUCTION 


A  pressing  problem  in  remote  sensing  is  the  detection  of  weak  targets  in  highly 
spatially  structured  optical  imagery.  Such  imagery  can  be‘generated  by  a  scanning  point 
detector  or  linear  photodiode  array,  or  by  a  staring  mosaic  photodiode  array.  Reviews  of 
current  image-processing  trends  have  illustrated  their  utility  for  enhancing  the  inherent 
information  content  found  in  remotely  sensed  multispectral  imagery  [l-b].  Most  of  these 
techniques  are  generally  suitable  for  feature  extraction  and  or  pattern-recognition 
applications,  and  as  a  consequence  are  not  always  helpful  in  the  detection  of  small,  weak 
targets.  Other  signal-processing  techniques  are  used  to  provide  target  detection  and 
localization  [7]. 

When  the  target  is  moving  with  respect  to  the  background,  several  signal-detection 
methods  can  be  applied  to  the  image  pixel  data  to  reduce  the  inherent  background  clutter 
in  the  imagery  and  enhance  target  detectability  [7-12].  For  example,  the  accepted  approach 
for  localizing  targets  in  staring  mosaic  data  sequences  is  to  frequently  bandpass  the  imagery 
by  using  either  an  analog  or  digital  filter;  e.g..  frame-to-frame  subtraction  or  a  third-order 
transversal  filter  [7 — 9j .  These  techniques  are  known  to  produce  excellent  results  provided 
the  target  is  relatively  strong.  Otherwise,  they  must  be  used  in  concert  with  other  signal¬ 
enhancing  methods  [7,  10-12], 

Mohanty  has  shown  that  potential  track  trajectories  of  a  weak  point  target  can  be 
assembled  through  exhaustive  search  of  all  possible  trajectories  in  a  given  frame  sequence 
|10].  Unfortunately,  for  a  staring  mosaic  array  of  reasonable  size,  the  number  of  possible 
trajectories  to  be  searched  can  be  large;  hence  it  is  an  unattractive  approach  to  target 
detection.  Barniv  used  dynamic  programming  to  reduce  this  approach  to  a  more  reasonable 
level  through  specific  knowledge  of  potential  speed  windows  and  target  shape  [l  l], 

A  more  robust  track  localization  technique  has  been  developed  by  Reed.  Gagliardi. 
and  Shao  for  detecting  weak,  moving  targets  in  background  clutter  by  using  three- 
dimensional  (3-D)  matched  filtering  [  1 2].  This  method  performs  moving-target-signature 
matched  filtering  in  the  Fourier  domain,  with  proper  signal  phasing  automatically  applied 
in  order  to  "coherently”  sum  the  target  energy  and  reduce  background  noise.  The  result  is  a 
set  of  matched-filter  peaks  indicating  detected  pieces  of  track.  Its  strengths  lie  in  its  ability 
to  operate  on  several  potential  target  tracks  simultaneously  (independent  of  their  respective 
space  time  origins)  and  to  bring  out  partial  as  well  as  lull  tracks  with  enhanced  signal-to- 
noise  ratios  often  substantially  in  excess  of  what  normally  would  be  obtained  from  spatially 
filtering  each  individual  frame  separately.  Its  weakness  is  that  this  method  requires  know¬ 
ledge  of  the  moving  target's  temporal  signature,  an  aspect  not  alwavs  known  a  priori.  This 
apparent  drawback  can  be  compensated  lor.  hut  at  the  cost  of  increasing  the  number  ot 
filter  banks  significantly  to  accommodate  all  possible  temporal  signature  variations  of  each 
prospective  target  of  interest. 

To  eliminate  this  last  cost  impact,  the  same  filtering  approach  can  be  reconfigured 
as  an  ensemble  moving-picture  processor  Its  specific  lorm  is  o!  a  space-shift  and  tinic- 
delay-and-integrate  matched  filter,  and  it  retains  much  of  the  processing  strength  described 
above.  However,  its  real  impact  is  to  reduce  the  number  of  filter  banks  to  the  level 
necessary  for  basic  target  velocity  discrimination,  and  this  can  be  reduced  further  bv  target 
parameter  estimation  by  means  of  data  fusion.* 


•Kendall.  W  B  .  Space  Computer  Corporation,  private  communication 


The  intent  of  this  report  is  to  further  examine  three-dimensional  matched  filtering 
and  to  show  how  it  can  be  applied  to  accomplish  target  detection  over  sequences  ol 
observed  scenes,  which  is  an  important  extension  of  classical  two-dimensional  image 
processing.  Examples  are  given  to  demonstrate  and  validate  these  results. 

2.0  ENSEMBLE  MOVING-PICTURE  PROCESSING 

Reed,  et  al.  1 2 J .  have  shown  that  the  optimum  three-dimensional  matched  filter  for 
a  moving  target  is  obtained  by  maximizing  the  received  signal-to-noisc  ratio  through  linear 
filtering.  Its  specific  form  is 


H(k:co) = 


A*(o)  +  k  •  v)  Sq( k ) 
k:co)  G( k ) 


where 


A(co)  =  jl  a(t)  exp  [- iuit  ;dt 

— OO 

+oo 

Sq( k )  =  J  sq(  r  1  exp  t-ik  •  r  }d~r 


<*>0(k:cL>)  =  j J  j  <Pq(£A)  exp  {- ik  •  r  -  icot> d “r « 


r,(k>  =  j  j  |li( _r  >|-  exp  -ik  *  _r  d~r 


where  the  authors  have  assumed  a  moving  target  of  spatial  profile  s()(r)  and  temporal 
intensity  function  a(t)  traversing  a  certain  portion  of  terrain  with  velocity  y.  This  terrain 
is  assumed  to  possess  colored  Gaussian  intensity  statistics  with  autocorrelation  0q( r:t ).  The 
entire  scenario  is  imaged  in  an  optical  system  with  point-spread  function  I h( r )| -  onto  a 
large  high-density,  high-speed  photodiode  array.  The  appendix  contains  a  detailed 
derivation  of  the  optimum  filter  H(k;co).  The  resulting  data  sequence  is  assumed  after 
sampling  to  be  quasi-stationary  in  intensity,  space,  and  time,  and  forms  a  continuous 
“moving  picture.”  which  will  be  image  processed  to  yield  the  existence  of  the  moving  target. 
On  the  assumption  that  the  signal  is  additive  to  the  noise,  the  two  possible  output 
sequences  are 


with  the  particular  choice  depending  on  whether  the  target  is  absent  or  present, 
respectively.  Here  nQ(r;t)  and  s(r:t)  are  the  received  background  and  target  intensity 
distributions,  respectively  (see  appendix  A). 

In  general,  the  exact  form  of  the  target  temporal  intensity  function  a(t)  is  not 
known  beforehand.  This  implies  that  the  size  of  the  matched-filter  bank  will  be  quite  large, 
since  it  must  encompass  all  possible  targets,  shapes,  velocity  windows,  track  orientations, 
and.  now.  intensity  functions.  Such  additional  complexity  would  generally  make  this 
approach  undesirable.  However,  if  the  above  three-dimensional  matched-filter  approach  is 
restructured  as  a  space-shift  and  time-delay-and-integrate  matched  filter,  significant  target 
signal-to-noise  ratio  gains  are  still  achievable,  and  with  reduced  processing  complexity.  In 
this  technique,  a(t)  is  replaced  as  the  temporal  intensity  variation  by  a  processing  window 
function  for  simultaneously  interacting  the  various  frame  groups  within  the  sequence.  The 
simplest  form  for  aft)  is  a  rectangular  time  function  of  width,  T  |  seconds,  centered  about 
t  =  0.  Mathematically,  this  is  expressed  as 


aft)  =  rect  [tl 


1  :  |t|<T,/2 
0  ;  otherwise 


As  will  be  shown,  the  input  sequence  must  be  greater  than  or  equal  to  T |  so  that  this 
technique  will  yield  optimal  results.  The  potential  signal-to-noise  ratio  in  this  case  is 
derived  below.  For  simplicity  and  without  loss  in  generality,  let  E[nQ(r:t)j  -  0.  Also,  let  the 
noise  spectrum  be  white  Gaussian  and  the  point-spread  function  of  the  optical  system 
| h(r)| 2  equal  5(r). 

For  an  image  sequence  of  long  duration,  the  input  signal  intensity  has  a  Fourier 
transform  of  the  form 

S(k  :co)  =  T0  sine  [<co  +  k  •  u)  Tq/2tt}  SQ(k  )  e  ~  ~  0  <  ' 


where 


Tq  =  frame  sequence  length  07  j )  =  MTq 
Tq  =  time  between  frames,  and 
sine  (7t.ct) 


sine  [a]  = 


ira 


(Si 


Figure  I  illustrates  a  target  moving  across  a  fixed  scene  in  M  =  Tq  Tq  frames  of  a 
sufficiently  large  sequence  of  images.  Using  equation  I  with  a  rectangular-window  signal 
function,  the  optimum  three-dimensional  filter  becomes 

Srf(  k ) 

H  (k:co)  =  Ti  sine  [Cco  +  k*  t>)Ti/2tr]  - —  i 

~  ~  N0 


for  a  white-Gaussian  background-clutter  target-detection  scenario.  This  implies  that  the 
output  image  has  a  mean  level  at  the  target  position  equal  to 


FINITE  DATA  SEQUENCE: 


Figure  1 .  Example  image  set  containing  a  moving  target  transvering  a  fixed  scene  during  a  finite  portion 
of  a  countably  infinite  data  sequence. 
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where 


Blip  { 


0  : 
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T1 
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From  equation  13.  it  is  apparent  that  the  matched-filter  pcak(s)  oeeur  within  the  eentral 
portion  of  the  second  image  sequence  of  time  length  1 0  Normally.  this  peak(s)  would  he 
expected  to  appear  towards  the  end  ot  the  data  sequence.  However,  the  eentral  processing 
peak  approach  was  chosen  out  of  a  desire  for  symmetry  in  the  processed  output  and  will  be 
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retained  for  the  remainder  of  this  report.  Both  output  approaches  are  correct,  but  those 
readers  desiring  the  classical  result  need  only  multiply  equation  9  by  the  appropriate  phase 
term  to  shift  the  filter  peak  to  the  sequence  end. 

The  image  variance  is  given  by 
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Thus  the  resulting  signal-to-noise  ratio  can  be  written  as 
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at  its  peak  location.  It  is  apparent  from  equation  20  that  three-dimensional  matched 
filtering  yields  a  peak  signal-to-noise  ratio  increase  over  two-dimensional  matched  filtering 
by  a  factor  of  T|.  the  temporal  width  of  the  rectangular  time  filter.  Clearly,  this  gain  can  be 
significant  for  most  reasonable  integration  times.  However,  it  should  be  emphasi'eH  T h * 
these  results  hold  only  for  constant  target  intensity.  Adjustments  in  the  required  SNR 
margin  should  be  made  for  intensity  \ariations  such  as  scintillation  and  lading  ol  the 
target,  which  might  occur  during  the  filter  integration  time.  In  the  following  section,  a 
sequence  of  computer-generated  images  with  a  modeled  airplane  mining  across  the  frames 
is  used  to  verifv  these  results. 


3.0  ENSEMBLE  MOVING-PICTURE  PROCESSING  OF  SIMULATED  DATA 


To  test  the  above  theoretical  formulation,  three  groups  of  30  sample  sets  were 
created.  Each  sample  set  consisted  of  19  white  Gaussian  images,  with  a  modeled  airplane 
propagating  across  each  of  them.  The  random  numbers  used  in  each  group  of  each  set  were 
independently  obtained  from  a  mathematical  procedure  suggested  by  Dillard  [|3J.  Each 
group  contained  zero  mean  image  sequences  with  an  average  noise  clutter  of  10.0.  The 
three  groups  are  distinguished  by  having  target  signatures  with  single-pixel  signal-to-noise 
ratios  of  6  dB.  0  dB.  and  -6  dB,  respectively.  Here  we  define  the  single-pixel  signal-to-noise 
ratio  to  be 


SNRj  = 


[S0-  B]- 


(21  ) 


with  Sq  —  intensity  value  of  the  target  within  one  pixel 

B=  mean  background  intensity  level 
=  0.0  for  this  simulation  i 


and  0=  Standard  deviation  of  the  background  clutter 
[=  10.0  for  this  simulation]. 


Figure  2  illustrates  an  example  frame  of  data  and  the  spatial  profile  of  the  simulation 
target.  In  this  figure,  the  spatial  dimensions  of  the  image  are  23  X  23  pixels  and  the  signal- 
to-noise  ratio  at  one  target  pixel  is  23.5  dB.  Figure  3  shows  the  same  frame  as  in  figure  2. 
except  that  the  target  signal-to-noise  ratio  is  0  dB.  The  speed  of  the  target  was  set  at  a 
uniform  I  pixel  per  frame.  Specific  matched  filters  were  generated  by  using  equation  I  for 
each  image  sequence  for  the  following  integration  time  windows:  T  =  I.  3.  4.  9.  II.  15. 
and  19.  Each  image  and  each  sequence  were  zero-filled  to  yield  45X45X45  arrays  to 
minimize  image  wraparound  when  processed  by  using  the  discrete,  finite-interval  Fourier 
transform.  Figures  4.  5.  and  6  plot  the  expected  signal-to-noise  ratio  given  in  equation  19 
for  the  matched-filter  peak  as  a  function  of  frame  number  for  initial  single-target-pixel 
signal-to-noise  ratios  (SNR)  equal  to  6  dB.  0  dB.  and  -6  dB.  respectively.  For  these 
calculations,  the  peak  signal-to-noise  ratio  can  be  written  as 

SNRp  =  10  log l0  T(  +  10  log10  ( 1 5)  +  SNRj  (dB)  (22a) 

=  10  log|0  IT,  )  +  11.76  dB  + SNR,  (22b) 

in  decibels.  The  first  term  represents  the  processing  gain  from  temporal  matched  filtering, 
the  second  represents  the  gain  from  spatial  matched  filtering  (equal  to  the  number  of 
independent  samples  composing  the  target),  and  the  third  term  is  the  inherent  signal-to- 
noise  ratio  ol  the  target  in  this  background  clutter.  These  last  two  terms  comprise  the 
expected  peak  signal-to-noise  ratio  obtainable  from  two-dimensional  matched  filtering. 
Using  equation  22b.  the  linear  increase  of  signal-to-noise  ratio  with  observation  time  can 
clearlv  be  seen  in  these  three  figures. 
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Figure  6.  Signal  to-noise  ratio  from  3  D  Mter.ng  SNR  -6  < j 8 


S' 

s 


I 


Figures  7  and  8  were  typical  filtered  images  obtained  by  three  matched  filters 
processing  the  image  sequences  from  which  figure  6  was  taken,  using  an  observation 
window  of  15  frames.  These  last  two  figures  were  created  by  squaring  the  real  and 
imaginary  parts  of  each  filtered  image  pixel  and  storing  their  individual  sum  as  the 
resulting  filtered  pixel  for  each  processed  image.  The  signal-to-noise  ratio  of  each  image 
was  calculated  by  sampling  the  highest  filtered  target  pixel  (the  peak  created  by  the 
matched  filtering),  calculating  the  first-order  moment  of  the  individual  frame  (minus  the 
target  area  contribution),  and  using  the  equation 


SNR 


M 


=  E{S^  +  N~)  _  j 
E{N2} 


(23) 


The  pixel  location  of  the  peak  intensity  was  in  row  14  and  in  a  column  position  equal  to 
the  frame  number  plus  two.  The  signal-to-noise  ratios  for  figure  7  were  around  23.5  dB. 
Specific  details  of  the  processed  results  are  as  follows. 

Figure  8  shows  resulting  signal-to-noise  ratios  for  the  filtered  6-dB  image  set  as  a 
function  of  the  observation  window  These  signal-to-noise  ratios  were  obtained  by  aver¬ 
aging  30  samples  of  measured  signal  levels  and  mean  noise  levels  within  a  I  X  I0X  I  X23 
window  of  each  image  and  using  equation  23.  It  is  apparent  from  this  figure  that  the  signal- 
to-noise  ratios  for  each  image  increases  with  integration  time,  with  the  greatest  gains  within 
the  center  images.  From  our  discussion  in  the  previous  section,  this  is  not  too  surprising 
However,  a  comparison  of  figure  8  with  figure  4  shows  that  the  measured  signal-to-noise 
ratios  have  a  more  graceful  falloff  from  their  peak  values  than  theory  predicts.  The  reason 
for  this  is  shown  in  figure  9.  Recall  that  figure  I  illustrates  tht  assumed  processing 
geometry  for  the  theoretical  development  given  above  We  have  an  infinite  (or  very  large) 
number  of  frames  in  which  a  certain  number  of  frames  will  have  a  target  propagating 
across  the  sequence  This  allows  maximum  buildup  of  signal  energy  and  a  constant  noise 
variance  to  occur  Figure  9  shows  the  real  situation  We  are  reallv  dealing  with  a  finite 
number  of  frames:  19  in  this  simulation  This  implies  that  the  noise  variance  will  be  subject 
to  the  same  energv  buildup  as  the  signal  energv  does,  and  will  be  given  bv 


VAR  |  vi  r  1 1)  =  Btti  j/ 
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This  is  because  the  ensemble  moving-picture  processing  is  a  sliding-interval  technique  and 
will  include  /ero-lilled  images  in  its  processing  chain  tor  large  integration  times  1  he  result 
is  a  filtered  signal-to-noise  ratio  equal  to 


0  dB  3-D  MATCHED  FILTER  RESULT 


Figure  9.  Realistic  image-processing  geometry  for  three-dimensional  matched  filtering. 
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rather  than  equation  19.  Figure  10  shows  a  comparison  of  equation  25  with  the  measured 
signal-to-noise  ratio  shown  in  figure  8.  It  is  apparent  that  these  results  agree  well  with  the 
revised  theory.  However,  it  is  also  clear  that  there  is  a  slight  oscillation,  on  the  order  of  a 
dB  or  less,  of  the  measured  signal-to-noise  ratio  for  the  smaller  filter  integration  times. 
Figures  1  I  and  12  compare  equations  12  and  24  with  the  measured  signal  and  noise 
variance  levels,  respectively.  It  is  clear  from  figure  I  I  that  the  measured  signal  level  follows 
the  theory  almost  perfectly,  while  figure  12  shows  an  oscillatory  nature  in  the  measured 
noise  variance  as  compared  to  theory.  This  may  be  a  result  of  our  using  an  infinite-interval 
processing  filter  with  finite-domain  Fourier-transformed  imagery,  whose  effects  average  out 
for  large  filter  integration  times.  However,  this  oscillatory  effect  is  small  and  does  not 
negate  the  power  of  the  technique. 

Figure  13  illustrates  the  measured  signal-to-noise  ratio  from  the  0-dB  initial  SNR 
simulation  with  equation  25.  As  noted  earlier,  the  image  sequences  used  were  drawn  from  a 
completely  different  random-number  sequence  than  was  used  to  create  the  6-dB  set.  It  is 
clear  that  these  results  also  agree  well  with  the  revised  theory,  with  the  same  oscillators 
behavior  at  the  lower  integration  times.  Figures  14  and  1 5  compare  equations  12  and  24 
with  the  measured  signal  and  noise  variance  levels,  respectively,  and  show  that  the 
oscillations  are  in  the  measured  noise  variance,  as  noted  before.  Figures  16.  I  7.  and  18 
repeat  the  same  comparisons  as  13.  14.  and  15.  respective!) .  except  lor  an  initial  signal-to- 
noise  ratio  of  -6  dB.  These  results  are  consistent  with  the  other  two  independent 
simulations  and  illustrate  that  three-dimensional  matched  filtering  can  bring  a  target's 
presence  out  ol  the  noise  in  a  predictable  fashion.  This  can  be  quite  important  lor  weak 
target  detection  in  strong  background  noise  or  when  target  lading  occurs 

Figures  ID.  13.  and  16  indicate  that  three-dimensional  matched  filtering  creates  a 
matched-filter  peak  within  one  frame  such  that  a  target  no-target  decision  can  be  made  h\ 
means  of  intensitv  thresholding  there  This  suggests  that  this  approach  mav  have  merit  in  a 
mov ing-target  indicator  (MTI)  scheme 


4.0  FNSEMBIF  MOVING-PIC Tl  RF  PROCESSING  AS 
A  MOVING-TARGET  INDIC  ATOR 


I  he  sensitivitv  of  ensemble  mm  ing-picture  processing  to  vclocitv  mismatch  can  be 
assessed  as  follows  From  our  development  in  section  2.  it  can  be  seen  that  a  vclocitv 
mismatch  will  onlv  affect  the  filtered  signal  level,  but  not  the  processed  noise  level  I  his 
implies  that  one  needs  onlv  to  recalculate  the  mean  signal  level  vtt  tl  I  he  input  sequence  '< 
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Figure  10.  Comparison  of  revised  theoretical  signat-to-noise  ratio  from  3-D  filtering  with  simulated 
data-processed  results.  SNR  =  6  dB. 
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Figure  12.  Comparison  of  filtered  image  noise  variance  with  revised  theory:  SNFT  =6  dB. 
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Figure  13.  Comparison  of  revised  theoretical  signal-to-noise  ratio  from  3-D  filtering  with  simulated 
data-processed  results:  SNR-  =0  dB. 
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Figure  14.  Comparison  of  measured  signal  level  with  procedural  results:  SNR:  =0  dB. 
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Figure  15.  Comparison  of  measured  noise-level  variance  with  revised  theory:  SNR(  =0  dB. 
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Figure  16.  Comparison  of  revised  theoretical  signal-to-noise  ratio  from  3-dB  filtering  with  simulated 
data-processed  results:  SNR-  =  -6  dB. 
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Figure  17  Comparison  of  measured  signal  energy  with  -evised  theory  SNR  -6  iB 
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In  the  above  equation'  v(|  i'  the  actual  target  veiocitv  v  |  i'  the  velocitv  'ettmg  ol  the  liltet 
and  bq  i'  the  Fourier  tran'lorm  ol  the  target  '  'patial  profile  in  it'  propel  angular 
orientation  !  he  mean  'ignal  level  alter  processing  i'  therefore 
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tor  v()  =  v  |  Equation  2*  'howv  that  the  re'ultrng  vignal  level  lor  a  velocitv  mismatch 
reduce'  the  inherent  'patial  matched-tiller  gain  hv  a  time-dependent  'patial  'hilt 
proportional  to  the  velocitv  dillerence  v  |  -  v()  I  he  sen'itiv  nv  ol  t h i '  mismatch  o  both  a 
function  ot  target  shape  as  well  as  velocitv  difference  hot  small  t.irget'  this  technique 
would  serve  as  an  VI  1 1  processor 

I  o  test  this  h v  pot hesis  the  b-d H  image  set  was  processed  with  a  set  ol  rn.it c  tied 
tillers  tuned  to  different  target  speeds  l  he  propagation  direction  was  assumed  to  be 
correct  \nv  directional  mismatches  would  reduce  the  tiller  signal  level  lurthet  so  oik 
might  consider  this  a  best  -  v  elocitv  -  mismatch  scenario  I  he  possible  target  speed'  were 
varied  from  0  t)  to  *  0  pixels  per  Irame  i  p  t  >  I  he  b-dH  target  traversing  the  input  tianic  set 
was  moving  at  a  constant  velocitv  of  I  pixel  pet  frame  Figure  1^  shows  the  measured 
signal-to-noise  ratio  at  the  SNR  peak  locations  used  to  create  figure  l<>  as  a  luuclion  of 
airplane  speed  1  he  tiller  observation  time  is  I1)  It  is  apparent  that  ensemble  processing 
with  incorrectlv  tuned  matched  tillers  vield'  lower  signal-to-noise  ratios  than  those 
resulting  Irom  the  use  ot  the  correct  tiller  as  expected  Interestinglv  there  is  a  notable 
asvmmctrv  between  the  processed  SNRs  depending  on  whether  tin  liltet  velocitv  o  ,ib,n. 
or  below  the  actual  target  velocitv  Speciticallv  the  resulting  signal-to-noise  tatios  degi.uk 
slightlv  more  slowlv  lor  those  (liters  whose  speeds  are  m  excess  ol  the  actual  speed  •>!  I  p  I 
This  is  more  clearlv  seen  in  figure  20.  which  depicts  the  SNR  for  the  peak  (tame  location, 
frame  10  From  equation  2K.  we  see  that  a  target  filter  mismatch  will  result  in  spatiai 
matched-filter  convolution  at  a  nonoptimum  offset  and  at  a  peak -time-conv  olm  u.n 
location  In  essence,  it  appears  that  the  resulting  filtered  target  sign. nine  is  essentia  I!  v 
optical  blurring  ol  the  expected  "matched"  fillet  peak  in  'he  t.uget  movement  loc.ito-n 
[  hat  is.  the  filter  can  coherently  stack  target  energy,  hut  in  an  ambiguous  \va\  I  his  cteateN 
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quite  desirable  from  a  data-processing  point  of  view  and  will  be  investigated  further.  If 
true,  this  means  that  3-D  filtering  for  mm ing-target  detection  could  be  viewed  as  a  coarse 
MTI  approach  and  could  hand  off  potential  target  information  to  other  systems  designed 
for  finer  target  localization.  Otherwise,  a  derivative  matched  filter  must  be  used  to  enhance 
shape  discrimination,  assuming  that  the  imagery  is  adequately  sampled  I4J  . 

5.0  CONCLUSIONS 

In  this  report,  we  have  shown  that  the  three-dimensional  filtering  techniques 
developed  by  Reed,  et  al.  12  ,  can  be  easily  modified  to  act  as  a  space-shift,  time-delay- 
and-integrate  matched-filtering  technique.  Computer  simulations  verified  the  validity  of 
this  modification.  The  results  demonstrate  signal-to-noise  ratio  gains  asymptotically 
growing  by  the  filter  integration  time  o\er  those  normally  obtained  by  two-dimensional 
matched  filtering  of  the  indisidual  frames  in  a  sequence.  The  specific  form  of  the  output  is 
ideally  suited  for  MTI  applications,  or  can  be  used  in  conjunction  with  some  other  target- 
speed  and  detection-localizing  algorithm. 
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APPENDIX:  DERIVATION  OF  OPTIMUM  3-D  MATCHED  FILTER 
FOR  COLORED  GAUSSIAN  NOISE 

Consider  the  remote  sensing  geometry  illustrated  in  figure  A-l.  An  optical  imaging 
system  is  located  a  distance  R  above  a  natural  terrain  with  temporal-spatial  intensity 
distribution  no'(r;t).  The  resultant  imagery  generated  by  this  sy  .cm  is  sampled  by  a  high- 
density.  high-speed  photodiode  array  and  stored  in  the  form  of  quasicontinuous  (both  in 
space  and  time)  image  sequences.  For  incoherent  background  illumination,  the  stored 
intensity  distribution  nQ(r;t)  can  be  shown  to  be  given  by 

+oo 

n0(r:t>  =  it  f  f  \h(r  -  t')\~  n'Q(r.t)  d~r‘  (A-l) 


where  x  is  a  real  constant  [l5]  .  In  this  equation.  h( r)  is  the  impulse  response  function  of  the 
optical  imaging  system  under  coherent  light  illumination.  The  magnitude  squared  of  this 
function  is  called  the  point-spread  function  (abbreviated  PSF).  and  is  linked  to  the  optical- 
transfer  function  (OTF)  of  the  imaging  system  through  the  relation 
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If  a  target  of  interest  moves  through  the  imaged  scene,  its  intensity  distribution  is 
also  mapped  onto  the  photodiode  array.  In  particular,  the  target  or  signal  intensity  has  the 
form 


s(  r;t)  =  x 
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where  s'  and  s  are  the  intensity  distributions  of  the  moving  target  before  and  after 
sampling,  respectively.  For  this  discussion,  let  us  assume  that  no  and  s  are  transferred  into 
storage  without  any  additional  noise  contamination.  This  implies  that  the  two  possible 
output  image  sequences  are 


y(_r;t)  =  nQ(^  ;t ) 


( A- 3 a  ) 


and 


y(  nt)  =  n0< jr ; t )  +  s(j;t) 
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Figure  A-1  Optical  remote  sensing  geometry. 


with  the  particular  form  ol  \ ( r:t )  depending  on  whether  the  target  is  absent  or  present, 
respectively  It  is  apparent  from  equation  A3-b  that  the  assumed  target  of  interest  is 
transparent  to  the  background  illumination,  is  blurred  to  the  point  that  this  additive 
relation  approximately  holds,  is  small  compared  to  the  pixel  area,  or  is  contained  in  frame 
to- Ira  me  ditlerence  data  Although  these  assumptions  may  not  include  all  possible  remote- 
sctiMng  scenarios,  the  application  ol  maximum-likelihood  detection  techniques  to  three- 
dimensional  image  processing  may  still  prove  useful  in  target-localuation  situations  not 
covered  bv  equation  A-3b. 


Given  the  above.  Reed,  et  al.  [12],  developed  the  optimum  three-dimensional  match 
filter  tor  moving-target  detection  shown  in  figure  A-2.  and  which  can  be  described  as 
follows:  The  target  of  interest  is  specified  to  possess  a  spatial  profile  sy(r)  and  a  temporal 
intensity  function  a(t).  In  addition,  it  traverses  the  imaged  terrain  at  a  velocity  v.  This  gives 
a  stored-target  intensity  distribution  of  the  form 


Figure  A-2.  3-D  matched-filter  processing  procedure. 
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with  rp  being  the  location  of  the  target  at  t  =  0.  Fourier  transforming  equation  A-4  yields 
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and 
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The  noise  in  the.received  pictures  is  given  hy 
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Fhis  implies  that  the  autocorrelation  function  tor  the  received  clutter  is  equal  to 
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The  power  spectrum  of  the  received  clutter  is  then 
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The  optimum  matched  filter  for  optical  target  detection  is  easily  shown  to  he 
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by  maximizing  the  received  signal-to-noise  ratio 
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